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Local  Approximation  on  Manifolds  Using 
Radial  Functions  and  Polynomials 

Jeremy  Levesley  and  David  L.  Ragozin 


Abstract.  The  main  focus  of  this  paper  is  to  give  error  estimates  for 
interpolation  on  compact  homogeneous  manifolds,  the  sphere  being  an  ex- 
ample of  such  a manifold.  The  notion  of  a radial  function  on  the  sphere  is 
generalised  to  that  of  a spherical  kernel  on  a compact  homogeneous  man- 
ifold. Reproducing  kernel  Hilbert  space  techniques  are  used  to  generate 
a pointwise  error  estimate  for  spherical  kernel  interpolation  using  a posi- 
tive definite  kernel.  By  exploiting  the  nice  scaling  properties  of  Lagrange 
polynomials  in  the  tangent  space,  the  error  estimate  is  bounded  above  by 
a power  of  the  point  separation,  recovering,  in  particular,  the  convergence 
rates  for  radial  approximation  on  spheres. 


§1.  Introduction 

There  is  currently  significant  interest  in  approximation  on  spheres,  related 
to  many  interesting  geophysical  problems.  There  are  a number  of  different 
approximation  methods  currently  available  on  spheres,  including  wavelets  [3], 
splines  [1],  and  the  subject  of  this  paper,  radial  functions  (sometimes  called 
spherical  splines)  [3,6].  Error  estimates  and  convergence  rates  for  radial  ap- 
proximation on  spheres,  of  an  optimal  nature,  are  recent  in  vintage  [5,4],  and 
rely  on  some  technically  demanding  mathematics.  In  this  paper  we  build  on 
an  idea  of  Bos  and  de  Marchi  [2]  in  order  to  provide  convergence  rates  for  ra- 
dial interpolation  on  a much  wider  class  of  manifolds:  the  reflexive,  compact 
homogeneous  spaces.  We  will  conclude  the  paper  by  proving  a local  spherical 
harmonic  polynomial  approximation  result  on  spheres. 

Let  Md  be  a d-dimensional  compact  manifold  with  a metric  d(-,  •)  which 
possesses  a transitive  group  G of  isometries.  The  group  acts  transitively  in 
that  for  every  x,y  £ Md,  there  exists  g £ G such  that  gx  = y.  If,  further- 
more, there  exists  g £ G such  that  gx  = y and  gy  = x,  then  Md  is  termed 
reflexive  (for  more  details  see  [11]).  Such  a manifold  is  a reflexive,  compact, 
metric,  homogeneous  space.  We  comment  that  we  can  always  embed  Md  in 
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some  higher  dimensional  Euclidean  space  Hd+r,  the  group  G being  a compact 
subgroup  of  the  isometries  of  !Rd+T'.  We  assume  that  the  metric  d(-,  •)  on  Md 
is  inherited  from  some  Euclidean  embedding. 

We  will  be  interested  in  interpolation  on  Md  using  continuous  zonal  ker- 
nels k( •,  •)  which  have  the  property  that  k{gx,gy)  = k(x,y ) for  all  x,y  g Md 
and  g € G.  Such  kernels  are  natural  generalisations  of  radial  functions, 
which  are  functions  only  of  distance,  which  is  itself  G-invariant.  Given  a 
set  { xi , ...  ,xn}  C Md  and  data  /i, . . . , /jv  6 IFt,  we  seek  a function  of  the 
form 

N 

sk{x)  = ^Taik(x,Xi) 

i=i 

such  that  Sk(xi)  = fi,  i = 1, . . . , N. 

Given  the  data  /;  = }{xf),  i = 1, . . . , N,  we  wish  to  bound  the  pointwise 
error  between  s/t  and  / at  x E Md.  We  make  no  assumption  on  the  data 
points  except  that  they  satisfy  a point  separation  criteria  in  some  subset  of 
Md  (see  Section  3). 

In  Section  2 we  will  introduce  some  necessary  harmonic  analysis  on  Md, 
discuss  the  notion  of  positive  definiteness  on  Md  in  brief,  and  give  a standard 
error  estimate,  which  we  will  use  in  Section  3 to  obtain  convergence  rates.  In 
Section  4 we  prove  a Whitney  type  error  estimate  for  local  spherical  harmonic 
approximation  on  the  sphere. 

§2.  Harmonic  Analysis  and  Error  Estimates 

For  a more  complete  version  of  the  brief  description  we  give  here,  see  [10,11]. 
Let  nd+r  be  the  degree  n polynomials  in  ftd+r,  the  space  in  which  Md  is 
homogeneously  embedded.  Then,  let  Tn  (the  spherical  polynomials  of  de- 
gree n)  be  the  restriction  of  these  polynomials  to  Md.  Furthermore,  let 
7in  '■=  'Pn  fl  'Pn-ii  where  the  orthogonality  is  with  respect  to  rf/t,  the  unique 
normalised  G-invariant  measure  on  Md: 

/ f9dli- 

J Md 

Then,  we  can  uniquely  decompose  Tin  into  irreducible  G-invariant  sub- 
spaces Enj,  each  of  dimension  dnj,  j = 1, ...  ,hn,  resulting  in  the  G-invariant 
decomposition 

L2(Md)  = ®-0®^15„J, 

Let  Xnj  be  the  orthogonal  projection  onto  Enj,  n = 0, 1, ... , and  j = 
1 ,...,hn,  and  7^nj(-,-)  be  the  kernel  of  this  projection.  We  will  consider 
interpolation  using  strictly  positive  definite  kernels  of  the  form 

OO  dn 

^ ^ 'y  ] anj'Rnj y)i 

n=0  j~\ 
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where  anj  > 0,  n = 0, 1, . . . , j = 1, . . . , hn,  and 

oo  dn 

EE  dnj  Q>nj  < oo.  (1) 

n= 0 j— 1 


We  will  approximate  functions  from  the  Hilbert  space 

oo  dn 

w = {fe  L2(Md)  : ll/ll2  :=  E < oo}, 

n=0  j=l 

where  ||  • ||2  denotes  the  L2{Md)  norm.  The  associated  inner  product  in  W is 


oo  dn 

(f,g)  :=  EE [Xnjf,Xnjg]/anj. 

n=0  j= 1 

The  condition  (1)  ensures  that  point  evaluation  is  a continuous  linear 
functional  in  W.  It  is  straightforward  to  show  that  k is  the  reproducing  kernel 
for  the  W:  f(x)  — (/,  k(x,  •)),/  € W,  x 6 Md . An  immediate  consequence  of 
the  reproducing  kernel  property  is  that  Sk  is  the  interpolant  of  minimum  W 
norm.  For,  if  g is  another  interpolant, 

N N 

( Sk  -g>sk)  = $>(/  ~ sk,Hxir))  = 5>(/(*0  “ = °* 

*=i  t=i 


Therefore, 

{9,9)  = {g-Sk  + sk,g-  Sk  + Sk)  = {9- sk,g- Sk ) + 2 (g  - Sk,sk ) + ( Sk,Sk ) 

= (g-sk,9-sk)  + (sk,sk),  (2) 

and  the  norm  minimisation  property  is  established.  Now,  following  the  stan- 
dard arguments,  see  e.g.  [8,9],  we  have,  using  the  fact  that  s*  interpolates  / 
at  xi, . . . ,xj\r, 


|/(x)  - sfc(x)|  = |(/  - sk,  k( x,  -)| 

N 

= |(/-s*,%,0  + E^%i’OI 

t = l 
N 

< ii/-^iii%o+E&%-oii 

*=i 

< 11/1111%,  o +!>%<,  oil 

t=i 
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for  arbitrary  /%,  i = 1, . . . , N,  where  we  have  used  (2)  in  the  final  step.  Our 
final  error  estimate  follows  from  the  fact  that 

N N N 

||fc(x,  •)  + Y^PMxi>-) II  = ikix, ')  + ^PMxir),  k(x , •)  + '}Tptk(xi,  -))5 

i= 1 t=l  t=l 

N N 

= {k{x,x)  -2^2(3ik(x,xi)  + 0iPjHxi>xi))^ ■ 

t = l ij  = l 


Defining 


N N 

P(x,Xi,...,xN)  :=  inf  (k(x,x)  - 2^  Pik{x,xi)  + Y]  PiPjKxi>xi))*  > 

pi.— iPnGK  — • . ^ 

t=i  t j=i 


we  sum  these  results  up  in 

Theorem  1.  Let  Sk  be  the  k-spline  interpolant  at  x\ , . . . , xn  6 Md  to  f £ W. 
Then,  for  every  x S Md, 

\f{x)  - Sfe(x)|  < ||/||P(x,  Xi,  . . . , in). 


§3.  Convergence  Rates  for  Radial  Kernels 

In  this  section  we  shall  give  pointwise  error  estimates  in  terms  of  the  point 
separation 

p :=  max  min  d{y,Xi), 
y£V  t=l,...,JV 

where  V C Md  contains  x,  the  point  at  which  we  are  measuring  the  error.  As 
we  shall  see  later  in  this  section,  producing  a pointwise  convergence  rate  from 
the  error  estimate  of  Theorem  1 requires  us  to  bound  Lagrange  polynomials 
related  to  a subset  of  the  interpolation  points.  Efforts  to  produce  convergence 
rates  on  the  sphere  foundered  because  it  is  difficult  to  bound  the  Lagrange 
polynomials  for  spherical  harmonic  interpolation  as  the  interpolation  set,  with 
a fixed  number  of  points,  scale  towards  x.  The  early  error  estimates  of  [3], 
of  0{p),  were  the  best  known  until  recently,  and  only  required  bounding  of 
the  constant  Lagrange  polynomial  for  a single  point.  Light  and  v.  Golitschek 
[4]  proved  boundedness  for  all  polynomials  on  Sd,  d > 2,  and  consequently 
achieved  0(pr)  approximation  for  radial  kernels  with  2r  continuous  derivatives 
on  the  sphere. 

A very  simple  proof  of  the  result  of  Light  and  v.  Golitschek  was  given  by 
Bos  and  de  Marchi  in  [2] . What  we  will  do  is  introduce  an  analytic  coordinate 
transformation,  and  construct  Lagrange  polynomials  in  the  tangent  space, 
which  is  a d-dimensional  Euclidean  space.  We  will  quote  a result  which  uses 
scaling  arguments  in  Euclidean  space  which  are  easy  to  perform,  observing 
that  distance  on  the  manifold  and  in  the  tangent  are  comparable. 
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Fig.  1.  Coordinate  chart  at  x. 

Let  x € V C M.  We  shall  assume  the  existence  of  a C°°-chart  (Uq,  ip) 
with  an  open  subset  U C Uo  satisfying  the  following  (see  Figure  1): 

1)  ip(U)  = V with  ip(0)  = x , 

2)  U0  = {y-z:y,z€U}. 

These  conditions  ensure  the  validity  of  Taylor  series  arguments  which  follow. 
Also,  since  U is  precompact,  if  is  bi-differentiable  and  the  metric  d is  assumed 
boundedly  equivalent  to  the  Euclidean  distance  on  lRd+’’, 

ci\\y  - z\\  < d{tp(y),ip(z))  < c2\\y  - z ||,  y,zeU. 

Let  v\ ,,vq  = V fl  X be  the  interpolation  points  in  V,  and  redefine 
p sup^gy  mini=it„.,Q  d(v,  vq).  Let  m = if~1(vi),  i — 1, . . . , Q.  Then,  from 
the  previous  equation  we  have 

rj  :=  sup  min  ||«  — Uj||  < p/c\.  (3) 

u£U  Q 

It  is  shown  in  [7]  that  provided  p and  hence  p are  sufficiently  small  to 
guarantee  that  Q > t,  we  can  make  a selection  of  interpolation  points  Vi, . . . ,vt 
(assuming  a convenient  ordering  of  the  points),  where  t = dim(ri2r_1),  such 
that  the  Lagrange  polynomials  pi,-..,pt  for  are  bounded  at  the 

origin: 

Pi(0)<Ct,  i = l,...,t,  (4) 

where  Cl  is  independent  of  i and  p.  Furthermore,  ui , . . . , ut  are  all  contained 
in  rjBb  :=  {rjw  : ||w||  < b}.  We  are  now  ready  to  prove  the  main  theorem  of 
this  paper: 
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Theorem  2.  Let  (j)  : H — > 3R  and  k(-,-)  = be  strictly  positive 

definite,  and  2 r-times  continuously  differentiable  in  each  variable.  Let  x £ 
S C M.  Suppose  that  the  interpolation  points  X\, . . . ,xn  satisfy 

sup  min  d(y,Xi)  = p. 

yeSi= l,-,n 

Then  for  all  sufficiently  small  p,  if  s*  is  the  k-spline  interpolant  to  f £ W, 
\f(x)-sk(x)\<C\\f\\pr, 
where  C is  independent  of  p. 

Proof:  First  we  now  choose  the  coefficients  /3i,...,pn  (appearing  in  the 
statement  of  Theorem  1)  as  follows.  Let  fi,  — pj{ 0)  if  xt  = ip(uj ) for  some 
j — 1, . ..  ,t.  Set  Pi  = 0 otherwise.  This  choice  of  coefficients  is  made  since 

t 

X)pj(°)9(u  j)  = q{°),  (5) 

j= i 

for  all  q G Then, 

t 

P(x,xi,...,xN)  < [^.(d(V>(0),V>(0)))  -2Y^Pj(0)<f>(d(t/’(0),i/>(uj))) 

t J=1  (6) 

«,J=1 

Since  <p(d(-,  •))  is  2r  times  continuously  differentiable  in  each  variable,  for 
fixed  w £ U we  may  expand 

= Y Ca(Z  - W)a  + R2riz)i  Z£U,  (7) 

|a|<2r 

where  R™r  is  a Taylor  series  remainder  satisfying 

Rl(z)<CR\\z-w\\2r,  (8) 

for  some  constant  Cr  independent  of  z and  w.  Putting  z = w in  the  above 
expansion  we  see  that 

<$*  = m-  (9) 

Putting  (7)  into  (6)  gives 

t 

P{x,XU...,XN)  < [<£(0)  — 2^Pi(0)^  Y C°a(U})a  + R2r(ujj) 

j= 1 |a|<2r 

+ Y Pi(°)Pj(°)(  U ca  («J  + («j))]  ' • 

i,j=1  |a|<2r 
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Using  the  polynomial  reproduction  (5),  we  get 

t 

P(x,x i,...,xN)  < j^-(?l)(0)-2^pJ(0)J?2rK) 

j-1 

+ EPi(°)(  E Co(-«i)“)  (10) 

i~  1 |a|<2r 

JL,  li/2 

+ E ft(°)pj(°)^(«j)J 

«>J=1 

where  we  have  used  (6)  and  (9)  in  the  above  argument. 

Now,  since  the  distance  function  is  symmetric,  for  any  u,w  € ]Rd, 

<f>(d(lp{u),  1p{w)))  = E Ca(u~  w)a  + R2r(u) 

\a\<2r 

= Y C“(®  - UT  + R2r  H- 
|a|<2r 

In  particular,  with  io  = 0we  get 

E £(-«)“  = E ci(ur+R°2r(u)-Ru o). 

|a|<2r  |a|<2r 

Substituting  the  last  equation  into  (10)  and  again  using  (5)  gives 

t 

p{x,x i,...,xN) < 

3= 1 

+ EM0)(  E C°(Ui)“+<(«)-^r(0)) 

*=1  |a|<2r 

1 !/2 

+ E Pi(0)P3(.°)R2r(uj)\ 
i,j= 1 

< [-E^(o)«K)  + ^(o)) 

3=1 

1 1/2 

+ E p^0^^0)^^)] 

i,j= 1 

t 

< CR{,  max  {|K -%-|U|Mi||})r  E IPi(°)ll2  + Pj(°)l 

t 

< CRTjr  Y CL( 2 + CL) 

i,j= 1 

<cP\ 
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using  (3),  (4)  and  (8),  recalling  that  u\,.  ,.,u t are  contained  in  an  rj  < p/ci 
scaled  ball  of  radius  b.  Substituting  this  result  in  Theorem  1 concludes  the 
proof.  □ 


§4.  Whitney-type  Estimates 

We  will  now  use  the  coordinate  system  suggested  by  Bos  and  de  Marchi  [2] 
to  prove  Whitney  type  estimates  for  approximation  using  spherical  harmonics 
on  Itd.  This  result  answers  a question  posed  by  L.  L.  Schumaker  during  the 
conference  for  which  these  are  the  proceedings.  For  x = (x2,  ■ ■ . , Xd+i)  € ]Rd, 
define  0 : ]Rd  ->  Sd  by  0(x2) . . . ,xd+i)  = ((1  - x?  )1/2,  x2,  ■ . . , xd+1). 
This  is  a smooth  parametrisation  of  a neighbourhood  ofei  = (1,0,...,  0)  € Sd. 

As  long  as  Y^= 2 X1  < s*n2  P'  *^en  d(ei,&(x))  < p. 

Without  loss  of  generality,  we  shall  consider  the  approximation  of  a func- 
tion / £ Ck(Sd),  using  spherical  polynomials,  in  any  spherical  neighbourhood 
Xp  of  e1;  where  maxy6xp  <L{e\ ,y)  = p < 7t/2.  In  fact,  we  will  prove 

Theorem  3.  Let  f £ Ck(Sd).  Then,  there  exists  a degree  k harmonic  poly- 
nomial pt-i  such  that,  for  every  ir/2  > p > 0, 

max \f(x)  -pfc_i(i)|  < C(f)pk, 
xexp 

where  the  constant  C(f)  does  not  depend  on  p. 

Proof:  The  crucial  element  of  this  proof  is  that  the  coordinate  mapping  0 
maps  polynomials  of  degree  k in  ]Rd,  the  tangent  plane  at  ei  coordinatised 
by  x2,...,xn+i,  to  polynomials  of  degree  k on  the  sphere.  Since  / o 0 g 
C'fc(0-1A: „),  we  can  perform  the  multivariate  Taylor  series  expansion 

Jo0(i)=  E cQxa  + Rk{f,x),  /n\ 

|a|<* 

where  Rk(f,x)  is  the  remainder  satisfying 

Rk(f,x)  < C{f)( max  |©_1(t/)|)fc.  (12) 

Letting  0 = 0(x),  and  defining  the  degree  k — 1 spherical  polynomial 

Pfc-iW  :=  E c«(0_1wr  = E 

|a|<fc  ja|<r 

equations  (11)  and  (12)  tell  us  that 

max  |/(0)  — Pk-i(0)\  < C(/)( max  |0_1(2/)|)fc. 

The  result  follows  because  maxyexp  |0_1  (?/)|  < p-  □ 
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